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METHOD AND APPARATUS FOR DETECTING ANOMALOUS SHADOWS 



BACKGROUND OF THE INVENTION 
Field of the invention 
5 The present invention relates in general to a method and 

apparatus for detecting anomalous shadows, and in particular 
to a method and apparatus for detecting anomalous shadows from 
a difference image, which is a comparative image representing 
the difference between two radiation i mages ■ 
m 10 Description of the Related Art 

K in a variety of fields, the comparative reading of two 

ffj or more images of the same subject obtained at different 

S photographing times to detect the difference between the 

^ aforementioned two or more images is performed, and an 

2 15 examination, etc, , of the subject of photographing is performed 

173 

p based on the detected difference, 

in the field of industrial product manufacturing, for 
example, an image obtained of a product in the new state is 
compared to an image of the same product that has been obtained 

20 after said product has been subjected to a durability test, and 
by focusing mainly on the area in which the largest difference 
between the two images appears, the portion thereof that should 
be improved to improve the durability of the product is 
determined. Further, in the field of medicine, by comparatively 

25 readi ng a pi u ral i ty of radi ati on i mages , whi ch have been obtai ned 
in a temporal series, of a diseased portion of a patient, a 
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physician can discern the course of the disease and the 
effectiveness of the treatment can be discerned and determine 
the course of future treatment. 

in this way, the comparative reading of two or more images 
5 obtained of the same subject is performed in a variety of fields 
on a daily basis, and in order to perform this comparative reading, 
there are cases for which these two or more images are outputted 
on an image display apparatus or the like. That is to say, the 
image is converted to luminosity and density signals, and then 
y3 10 di spl ayed on an i mage di spl ay apparatus or the 1 i ke , or outputted 
ID by a printer to a film or other such medium. 

CP For cases in which the two or more images that are the 

S3 objects of comparative reading are to be displayed, etc., 

although the general practice is simply to display said images 
ni 15 next to each other, the area of the images of most interest to 
O an operator performing a comparative reading thereof is the area 

showing the difference between said images. However, in a case 
in which, for example, two images that are to become the objects 
of comparative reading as described above are simply displayed 
20 next to each other, as the size of the difference between becomes 
smaller the detection thereof becomes more difficult; therefore, 
there is a demand for an improvement in the performance of 
comparative reading. 

Therefore, interimage computational processes, starting 
25 with a subtraction process between the pixels corresponding to 
each of the respective two or more images that are the objects 
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of comparative reading, is performed, and the difference 
therebetween is enhanced. Because it becomes possible to 
accurately render the difference between the images visually 
discernable, by enhancing only the difference between the images 
5 in this way, to an operator performing the comparative reading, 
it is possible to prevent the overlooking of the difference 
between the images, particularly in the medical field with 
respect to the course of a di sease or a di seased portion requi ri ng 
treatment. For example, on a difference image (a so-called 
5" 10 temporal subtraction image) obtained by the above-described 
JS{ interimage processes and representing the difference between 

% two chest X-ray images obtained in a temporal series, it is 

£ possible to display even small, early-stage lung tumors, which 

are substantially 100% treatable, and the extraordinary 
ns 15 effectiveness for the use in diagnosis is widely acknowledged. 
% Further, a method wherein, by adding the above-described 

i== subtraction image to a temporal image, the change that has 

occurred in a diseased portion becomes easier to discern visually , 
has also been proposed (refer to Japanese unexamined Patent 
20 Publication No. 8(1996)-77329) . 

However, although the process of selecting the diseased 
portion on the difference image is performed visually and 
manually by an operator performing a comparative reading as 
required, the accuracy of said selection is dependent upon the 
25 experience level and ability of the operator, and it is not a 
foregone conclusion that said selection will always be an 
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objectively viable selection. In a case, for example, in which 
the objective of radiation images that are to be the objects 
of comparative reading is the detection of lung tumors, it is 
necessary to extract the small, faint suspected anomalous 
5 shadows that are one indicator of cancerous growths; however, 
it is not guaranteed that those shadows will be accurately 
sel ected . Therefore , there i s a demand for a method and apparatus 
capable of accurately detecting the anomalous shadows that is 
not dependent upon the skill level of the operator. 

10 In order to respond to this demand, CADM (Computer Aided 

Diagnosis of Medical images), which has an objective of 
accurately and automatically detecting the suspected anomalous 
shadows by employing computational processes, has been 
advancing in recent years. 

15 That is to say, CADM technology, by performing the 

detection of the above-described suspected anomalous shadows 
based on performing a detection process for detecting the density 
distribution and shape characteristics thereof utilizing a 
computer, is a technology that automatically detects the 

20 suspected anomalous shadows. A variety of CADM technologies for 
application to chest X-ray images and which provide for a high 
probability of detecting diseased portions such as lung tumors 
have been proposed (refer to, for example, Japanese Unexamined 
Patent Publication No. 9(1997)-185714) . 

25 Here, the method proposed in the aforementioned Japanese 

Unexamined Patent Publication No. 9(1997)-185714 is a method 



4 



for applying CADM technology to an original image and an energy 
subtraction image (that is , the difference image between a signal 
enhanced image and a signal suppressed image) for detecting 
suspected anomalous shadows. However, there is a possibility 
5 that the early-stage 1 ung tumors in an original image or an energy 
subtraction image will be overlooked even by an experienced 
diagnostician; even by applying CADM technology, it has been 
nearly impossible to extract the suspected anomalous shadows 
of diseased portions such as early-stage lung tumors, which are 
O 10 easily overlooked- 

£ Therefore, applying CADM technology to the 

2 above-described difference image has been considered. Here, 

^ when the above-described difference image is to be formed, 

: although it is necessary to match the positions of two images, 

J 3 ! 15 even if position matching is performed, slight misalignments 
b? occur at the outlines of the structuring elements (positional 

^ misalignments), and artifacts due to these positional 

misalignments appear in the difference image. In particular, 
if the subject of photographing is a living-tissue subject, the 
20 3-dimensional misalignments due to fluctuations in the angle 
between the irradiation direction of the radiation and the body 
axis of the subject of photographing, the angle of entry of the 
radiation relative to the orientation of the subject of 
photographing, etc. , that appear in a 2-dimensionally displayed 
25 image are exceedingly difficult to correct (position matching) 
therefore, artifacts remain in the difference image. 
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Essentially, because only the actual difference portion 
representing the change, etc. of the diseased portion appears 
if artifacts do not occur in a difference image, by detecting 
the image portion that appears in a difference image, the size, 

5 etc. of the position of the diseased portion can be efficiently 
determined. However, if artifacts appear in the difference image, 
the above-described actual difference between the images does 
not stand out within the difference image due to the presence 
of the artifacts, and even if CADM technology is applied, there 

10 is a fear that it will not be possible to detect the difference. 

SUMMARY OF THE INVENTION 

The present invention has been developed in view of the 
circumstances described above, and it is a primary object of 
the present invention to provide an anomalous shadow detecting 

15 method and apparatus capable of adequately detecting anomalous 
shadows within a difference image, even for cases in which there 
are artifacts present in the within the difference image. 

According to the suspected anomalous shadow detecting 
method and apparatus of the present invention, by subjecting 

20 the difference image to an image process wherein the difference 
between the two images on which the difference image is based 
is enhanced relative to the artifacts, and the artifacts, which 
can have a negative effect on the diagnostic reading of the 
difference image, are reduced relative to the actual difference 

25 between the aforementioned two images, which is the area of 
interest in the image. 
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That i s to say , the method of detecti ng suspected anomal ous 
shadows according to the present invention comprises the steps 
of: obtaining a difference image representing the difference 
between two images, which have been obtained of the same subject 

5 at different photographing times, by subjecting said two images 
to an interimage process for obtaining the difference between 
said two images; obtaining a processed difference image by 
subjecting the obtained difference image to an image process 
wherein the actual difference between the two images is enhanced 

10 relative to the artifacts appearing due to misalignment of the 
positions of the structural elements thereof; and detecting the 
actual difference between the aforementioned two images from 
the processed difference image as suspected anomalous shadows. 
Here, as to the aforementioned two or more images, 

15 radiation images that have been obtained of the same subject 
in a temporal series, each at a different time, and which become 
the objects of a comparison that is performed to determine the 
temporal change can be applied, and the present invention is 
particularly effective when the subject of photographing is a 

20 living body. Because each of the internal structural elements 
of a living body move in a 3-dimensional manner, the relative 
positions of the structural elements between two images is easily 
thrown off due to the occurrence of positional misalignments. 
As to the interimage process, it is preferable that a 

25 subtraction process corresponding to the structural positions 
between the two images is applied. This is due to the fact that 
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by representing the difference between the two images as the 
difference between image signals, it is simple to render the 
difference clearly. However, the interimage process is not 
limited to being a subtraction process; the application of a 

5 division process or the like corresponding to the positions of 
the structuring elements within two images cannot be ruled out. 
As to the subtraction process, a simple subtraction computation 
or a weighted subtraction computation can be employed. An 
interimage calculated image obtained by a subtraction process 

10 is generally referred to as a subtraction image. These 

subtraction images include: an energy subtraction image (simple 
subtraction or weighted subtraction) obtained based on two 
original images having mutually different energy distributions 
(the two original images: a high energy exposure image <a normal 

15 radiation-image>, and a low energy exposure image <a 

high-pressure suppressed image>) and which have been obtained 
at substantially the same time in terms of a temporal series; 
a temporal subtraction image obtained based on two original 
images obtained in a temporal series, each obtained at different 

20 times; and a DSA (Digital Subtraction Angiography) obtained of 
images photographed of a blood vessel both before and after the 
injection of a dye; etc. 

As to the image process for enhancing the actual difference 
between two images relative to the artifacts, an image process 

25 that enhances the actual difference more than the artifacts, 
or an image process that suppresses the artifacts more than the 
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actual difference can be applied. 

As to the image process for suppressing the artifacts 
relative to the actual difference between two images, it is 
preferable that a morphology process based on a morphology 
5 computation employing structuring elements that are larger than 
the artifacts or smaller than the actual difference is applied. 
This is due to the fact that the artifacts remaining in a 
difference image can be effectively suppressed with respect to 
the aforementioned actual difference. Note that the "morphology 
10 process" in general , has been developed as a set theory occurring 
2j in an N-dimensional space, and is often applied to images, which 

5 are 2-dimensional spaces (refer to Japanese Unexamined Patent 

j| Publications. 8(1996)-272961, 9(1997)-248291, 9(1997)-91421, 

u etc.). Here, a simple explanation of this morphology process 

jij 15 will be given below, using a concentration image as an example. 
q A concentration image is regarded as a 3-dimensional space 

constituted of a certain number of points (x, y) having 
respective heights corresponding to a concentration value f (x, 
y) ; a 1-dimensional function f (x) can be considered appropriate 
20 for this area. As shown in formula (1), the structuring element 
g utilized in the morphology process is a symmetric function 
that is symmetrical about an origin point. 

g s (x)=g(-x) (1) 
The value is 0 at the defining area, and this defining area G 
25 is defined as shown in the following formula (2). 
G={-m, -m+1,..., -1,0,1,. ..,m, -m} (2) 

9 
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At this time, the basic form of the morphology process 
is as shown formulas (3) - (6), and is an extraordinarily simple 
computation . 

dilation; [f®G 0 ](i)=max{f (i-m),.-^ (i),-,f (i + m)} (3) 

erosion ;[fQG 9 ](i)=min{f (i-m),-,f (i),-,f (i + m)} (4) 

opening ;f g -(f0g 0 )®g (5) 

closing ;f a = (feg 0 )eg (6) 

That is to say, the dilation process is an operation that 
searches for the maximum value within the range of the width 
±m (a value defined by the structuring element and equal to the 
mask size shown in Fig -19a), which has the pixel of interest 
as the center pixel thereof, (refer to Fig. 19A) . On the other 
hand, the erosion process is an operation that searches for the 
minimum value within the range of the width ±m, which has the 
pixel of interest as the center pixel thereof (refer to Fig. 
19b). Further, the opening-process is an operation consisting 
of performing the erosion process followed by the dilation 
process; that is to say, an operation that searches for the 
minimum value first and then searches for the maximum value. 
On the other hand, the closing- process is an operation consisting 
of performing the dilation process followed by the erosion 
process; that is to say, an operation that searches for the 
maximum value first and then searches for the minimum value. 

in other words, the opening-process smoothes the 
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density-curve f (x) from the low-density side thereof, and is 
appropriate for controlling the up-pointing peaks of density 
fluctuation portions that fluctuate within the range spatially 
narrower than the mask, which has a size of 2m (refer to Fig. 
5 19C) . 

On the other hand, the closing-process smoothes the 
density-curve f (x) from the high-density side thereof, and is 
appropriate for controlling the down-pointing peaks of density 
fluctuation portions that fluctuate within the range spatially 
5 10 narrower than the mask, which has a size of 2m (refer to Fig. 19d) . 
jj Here for cases in which a signal becomes a high-density 

S high-signal level as the density value thereof increases, 

S| because the size relation with respect to the case in which the 

y, i mage si gnal val ue of the densi ty val ue f (x) i s a hi gh-1 umi nosi ty 

JJ 15 high-signal level is reversed, the dilation process for a 
54 high-density high-signal level signal and the erosion process 

for a high-luminosity high-signal level signal (refer to Fig. 
19B) are identical, and the erosion process for a high-density 
high-signal level signal and the dilation process for a 
20 high-luminosity high-signal level signal (refer to Fig. 19A) 
are identical; the opening process for a high-density 
high-signal level signal and the closing process for a 
high-luminosity high-signal level signal (refer to Fig. 19d) 
are identical, and the closing process for a high-density 
25 high-signal level signal and the opening process for a 

high-luminosity high-signal level signal (refer to Fig. 19C) 

11 



are identical . 

Therefore, by subjecting an image signal representing an 
original image to an opening or closing morphology process in 
this manner, granularity 0" .e. , the noise occurring as an image 

5 signal) can be suppressed (or eliminated) (refer to, for example, 
"Morphology", by obata, Corona Press). Note that for the sake 
of simplicity an explanation of a case in which a 1-dimensional 
mask (a structural element) has been employed in morphology 
processing of a 1-dimensional density distribution has been 

10 given, however, in applying a morphology process to an image 
extending 2-dimensionally , a plurality of 1-dimensional 
structuring elements can be set within the surface of the 
2-dimensional image and the process can be performed a plurality 
of times, that is, once for each respective structural element, 

15 or a 2-dimensional structural element can be set, and the process 
can be performed once applying this 2-dimensional structural 
element. 

Note that according to the anomalous shadow detecting 
method of the present invention, it is preferable that the 
20 substantially round-shaped differences from among the actual 
differences appearing in the processed difference image are 
detected as the suspected anomalous shadows. 

The anomalous shadow detecting apparatus according to the 
present invention is an apparatus for implementing the anomalous 
25 shadow detecti ng method of the present i nventi on , and compri ses : 
an interimage processing means for obtaining a difference 
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image representing the difference between two images which have 
been obtained of the same subject at different photographing 
times, by subjecting said two images to an interimage process 
to obtain the difference between said two images; 

5 an image processing means for obtaining a processed 

difference image by subjecting the aforementioned difference 
image to an image process wherein the actual difference between 
the two images on which the aforementioned difference image is 
based is enhanced relative to the artifacts appearing due to 

10 misalignment of the positions of the structural elements 
thereof; and 

a detecting means for detecting the actual difference 
between the aforementioned two images from the processed 
difference image as suspected anomalous shadows. 

15 as to the aforementioned two or more images on which the 

difference image is based, implementing the present invention 
is more effective and advantageous if applied to radiation images 
that have been obtained of the same subject in a temporal series, 
each image being obtained at a different time, and which are 

20 to become the objects of compari son for determi ni ng the temporal 
change therebetween. 

As to the interimage process, it is preferable that a 
subtraction process corresponding to the structural positions 
between the two images is applied. 

25 as to the image processing means, a means for carrying 

out a process which suppresses the artifacts relative to the 

13 
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actual difference between two images, or conversely, a means 
for carrying out a process which enhances the actual difference 
between two images relative to the artifacts can be applied. 
For cases in which the image processing means is a means for 
carrying out a process which suppresses the artifacts relative 
to the actual difference between two images, it is preferable 
that a morphology process based on a morphology computation 
employing structuring elements that are larger than the 
artifacts or smaller than the actual difference is applied. 

Note that as to the detecting means, it is preferable that 
a means for detecting the substantially round-shaped 
differences from among the actual differences appearing in the 
processed difference image as the suspected anomalous shadows 
is employed. 

According to the method and apparatus for detecting 
anomalous shadow according to the present invention, by 
subjecting an interimage image to an image process wherein the 
actual difference between two images is enhanced relative to 
the artifacts appearing in the interimage image due to 
misalignment of the structural positions between said two images, 
the artifacts, which have a negative effect with regard to the 
detection of the suspected anomalous shadows can be reduced 
relative to the actual differences between the two images, which 
are regarded as the suspected anomalous shadows. Accordingly, 
the suspected anomalous shadows can be detected accurately, with 
no negative effect from the artifacts. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is an overview drawing of the configuration of 

a network including an embodiment of the suspected anomalous 

shadow detecting apparatus according to the present invention, 
Figures 2A and 2B illustrate two original images that are 

to be subjected to a temporal subtraction process by the 

interimage processing apparatus, 

Figures 3A and 3b are drawings illustrating global 

matchi ng , 

Figures 4A and 4B are drawings illustrating local 
matchi ng , 

Figure 5 is a drawing illustrating a non-linear transform, 
Figures 6A, 6B, and 6c are drawings illustrating a first 
original image, a second original image and a subtraction image 
for a case in which the positions of the first and second images 
have been completely matched, 

Figure 7 is a drawing illustrating an example of an actual 

subtraction image, 

Figure 8 is a drawing illustrating a temporal series 
according to the present invention, 

Figure 9 is a drawing illustrating the division of a 
subtraction image into small regions, 

Figure 10 is a drawing illustrating the image within a 
selected small region, 

Figure 11 is a drawing illustrating the image within a 
selected small region in the state in which the center has been 
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set, 

Figure 12 is a drawing illustrating a coordinate image 

of selected small region, 

Figure 13 is adrawing illustrating an edge detection mask, 
Figure 14 is a drawing illustrating an edge image, 
Figure 15 is a drawing illustrating the histogram of an 

edge portion occurring in an edge image, 

Figures 16A and 16b are drawings illustrating template 

matching, 

Figure 17 is a drawing illustrating the detection state 
of a small region, 

Fi gures 18A and 18b are a drawi ngs ill ust rati ng the di spl ay 
state of a subtraction image in which the suspected anomalous 
shadows have been detected, and 

Figures 19A, 19b, 19C, and 19d are graphs each illustrating 
the concept of a morphology process. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Hereinafter, a preferred embodiment of the method and 
apparatus for detecting anomalous shadows according to the 
present invention will be explained with reference to the 
attached drawings. 

Fig. I is a drawing of a medical image network 100 including 
an anomalous shadow detecti ng apparatus, which is an embodiment 
of the method and apparatus for detecting anomalous shadows 
according to the present invention. 

connected to the network 100 shown i n Fi g . 1 i s an anomal ous 
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shadow detecting apparatus provided with, for example: a CT 
(Computer Tomography imaging) apparatus; an mri (Magnetic 
Resonance Imaging) apparatus; a CR (Computer Radiography) 
apparatus 50 or other medical image forming apparatus; a data 
5 base 70 for cumulatively storing each type of diagnostic medical 
image formed by these medical image forming apparatuses; an image 
display apparatus 20 for displaying images temporarily stored 
in the data base 70 or images that have been transmitted directly 
from the image forming apparatuses; an interimage processing 

10 means 11 for performing a temporal subtraction process, which 
is one example of an interimage process, based on two or more 
image data input thereto and obtaining a temporal subtraction 
image (an example of an interimage processed image, hereinafter 
referred to as a subtraction image Psu) ; and an image processing 

15 means 12 for subjecting the subtraction image Psu formed by the 
interimage processing means 11 to the process to be described 
below and; a detecting means 13 for detecting the suspected 
anomalous shadows from the temporal subtraction image that has 
been subjected to the image process (hereinafter referred to 

20 as a processed subtraction image Psu"). Note that a variety of 
other image outputting devices (not shown) such as a printer 
for outputting images to film or the like are connected to the 
network. 

The CR apparatus 50 records the radiation image of a 
25 subject of photographing on a stimulable phosphor sheet having 
a layer containing stimulable phosphors by irradiating said 
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stimulable phosphor sheet with the radiation energy that passes 
through said subject of photographing upon the irradiation 
thereof with a radiation, and then irradiates said stimulable 
phosphor sheet with a laser light or other excitation light and 
5 photoelectrical^ reads out the stimulated emission emitted 
thereupon as a quantity of light corresponding to the radiation 
energy that has been recorded on the phosphor sheet, whereby 
a radiation image of the radiation energy that has passed through 
the subject of photographing is obtained as a digital image; 

10 these CR apparatuses are widely used in hospitals and other 
medi cal f aci 1 i ti es - 

The QA-WS (a workstation for image quality control) 60 
disposed between the CR apparatus 50 and the network 100 checks 
the quality of an image formed by the above-described CR 

15 apparatus 50 or other image forming means, and is a work station 
provided with a function for outputting a command to an image 
forming apparatus (the CR apparatus 50, etc.) to re-obtain an 
image, etc., as required- The QA-ws 60 employed in the current 
embodiment is a means provided for displaying an image P formed 

20 by the CR apparatus 50 before storing said image Pin the data 
base 70, and checking the image qualities such as the image 
density, the contrast, etc., as well as the photographing 
envi ronment . 

The image display apparatus 20 is not simply a means for 
25 displaying as a visible image an image inputted thereto via the 
network 100; said image display apparatus 20 is also provided 
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with a function to serve as an input terminal for inputting a 
command to perform a subtraction process on two or more images 
P obtained of the same subject in a temporal series, etc. 

The interimage processing means 11, as described above, 
5 forms a subtraction image Psu, based on two or more chest X-ray 
images Pi and P2, each image having been obtained of the same 
subject at a different point in time in a temporal series and 
which have been inputted thereto from the database 70, which 
represents the difference between these two images Pi and P2; 

10 the forming of this subtraction image Psu is carried out by fi rst 
globally matching between the images Pi and P2 the structural 
elements (the sternum, the thorax, the lungs, the heart and other 
organs, etc.) appearing within each of images Pi and P2, and 
then dividing the images to obtain a plurality of small local 

15 regions and performing local matching so as to correlate the 
positions of each of said local regions. Then a subtraction 
process corresponding to the thus correlated pixels is performed 
between the two images obtained by this 2-step matching process. 
Note that the particulars of each position matching process are 

20 described below. 

The content of the image process that the image processing 
means 12 performs on the subtraction image Psu formed by the 
interimage processing means 11 consists of subjecting the 
artifacts L, which are caused by misalignments between the 

25 correlative positions of the structuring elements of the image 
Pi and P2, appearing within the subtraction image Psu, which 
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represents the difference between the images Pi and P2, to an 
image process that suppresses said artifacts L relative to the 
actual difference K, which is the portion between the two images 
Pi and P2 in which the density and the outline of the forms differ 
5 not because of a misalignment between positions, but due to the 
presence of a shadow of a tumor or diseased portion in only one 
of the two images Pi and P2 , or due to the temporal change in 
a diseased portion due to the treatment thereof or the 
progression thereof (herein, after referred to simply as a 
O 10 diseased portion) appearing in both images. 

;D More specifically, because there are numerous cases in 

5( which the artifacts L appearing in a subtraction image Psu are 

J£| characterized in that said artifacts L have an elongated shape 

I, and are located along the outlines of the structural elements, 

tfj 15 although the diseased portion K, which is the actual difference 
;2: between the two images Pi and P2, fluctuates in size 

H? corresponding to the course of the disease, the diseased 

portion(s) K can be recognized as the region(s) of a shape having 
a width which expands into a more substantially round or 
20 rectangular shape, etc., than that of the elongated shape of 
the artifacts L. Accordingly, by applying a morphology process 
utilizing a structuring element that is smaller than the diseased 
portion K and has a width(the length in the direction 
substantially perpendicular to the lengthwise direction of the 
25 artifacts L) larger than that of the artifacts L (such as a round 
structuring element having a diameter of said size, or a 
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rectangular structuring element having a side of said size, etc.) 
as the aforementioned image process to which the subtraction 
image psu is to be subjected, a diseased portion K of a size 
larger than that of the structuring element can be caused to 
5 remai n , whi 1 e the arti facts L of a 1 ength shorter than the 1 ength 
thereof in any direction can be eliminated. 

Note that for cases in which a morphology process is 
applied to a subtraction image Psu, which is a concentration 
image, the subtraction image Psu is separated into a high-density 

10 image Psul formed of the portions of the image having a density 
equal to or higher than the base density (the density of the 
portions of the images Pi and P2 of which the difference between 
the respective densities thereof is 0) (the entire image Psul 
is formed of said image portions having a density equal to or 

15 greater than the base density, and the low-density portions, 
which have been converted to the base density) , and a low-density 
image Psu2 formed of the portions of the subtraction image Psu 
having a density equal to or less than the base density (the 
entire image PSU2 is formed of said image portions having a 

20 density equal to or less than the base density, and the 

high-density portions, which have been converted to the base 
densi ty) ; each of the hi gh-densi ty i mage Psul and the 1 ow-densi ty 
image Psu2 are subjected to the above-described morphology 
process, and because it is preferable that the two images 

25 obtained thereby are subjected to an addition process (in which 
the images are weighted at a ratio of 0.5 to 0.5), hereinafter, 
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the current embodiment will be explained in terms of said 
addition process. 

The detecting means 13 detects the diseased portion K from 
the processed subtraction image Psu', in which the artifacts 

5 L have been suppressed, as an anomalous shadow(s) . More 
specifically, the diseased portion(s) K is detected as an 
anomalous shadow(s) by use of a template matching method 
utilizing a template having the shape of the diseased portion 
K (substantially round or substantially rectangular): refer to 

10 Japanese unexamined Patent Publication No. 9(1997) -185714, in 
which a method for detecting the diseased portion K as the 
anomalous shadow(s) by use of a neural network is disclosed. 

Next, the operation of the suspected anomalous shadow 
detecting apparatus according to the current embodiment will 

15 be explained. 

A series of chest X-ray images Pi, P2, Pn, each image 
having been obtained at a different point in time, of a specific 
patient that have been obtained in advance by the CR apparatus 
50 is subjected to an image quality check by the QA-WS 60, and 

20 then output to the data base 70 and cumulatively stored therein. 
First, an example in which a command to display a 
subtraction image Psu based on an image Pi (the first image), 
which is the image within a temporal series of images obtained 
of the same subject that has been obtained at the earliest point 

25 in time of said temporal series, and an image P2 (the second 
image), which is an image that has been obtained at a point in 
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time later than that of said image Pi, has been issued via the 
image display apparatus 20, which also serves as an input 
terminal for inputting commands, will be described. 

The command inputted to the image display apparatus 20 

5 is inputted to the data base 70 via the network 100. The data 
base 70 inputs, according to the command received from the image 
display means 20, the first image Pi and the second image P2 
to the interimage processing means 11 via the network 100. 

The interimage processing means 11 first performs global 

10 matching between the positions of the structural elements 

appearing in the images Pi and P2, and then divides said images 
Pi and P2 into a plurality of small local regions and performs 
local matching so as to correlate the positions of each of said 
local regions. That is to say, taking the image P2 as the base 

15 image from between the images Pi and p2 shown in Figs. 2A and 
2B, respectively, the first image Pi is subjected to an affine 
transform, and is globally matched to the newer image from among 
the temporal series, P2 (the second image), (refer to Figs. 3A, 
3B). continuing, the second image P2 from among both the images 

20 Pi and P2 that have been globally matched is subjected to a 
process employing a template and is divided thereby into a 
plurality of virtual small rectangular template regions T2. Then, 
the first image Pi is correlated to each of the template regions 
T2 of the second image P2 , a search region Rl larger than each 

25 template region T2 is set, and for each search region Rl, a 
corresponding and substantially matching template region T2 of 
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the second image P2 is obtained (refer to Figs. 4A and 4B) ; based 
on the corresponding positional relation between each template 
region T2 of the second image P2 and each corresponding template 
region of the first image Pi, after each corresponding template 
5 region of the first image Pi has been globally matched to each 
respective template region T2 of the second image P2 so as to 
be substantially matched, the enti re fi rst image Pi is subjected 
as a whole to a non-linear transform (refer to Figs. 5A and 5B) , 
and the positions of the two images Pi and P2 are again matched. 

10 Further, the interimage processing means 11 correlates 

the pixels of the base image P2 and the first image Pi that has 
been subjected to the transform of the above-described 2-step 
position matching process, and subtracts said first image Pi 
from said second image P2 to obtain a subtraction image Psu. 

15 Here, for a case in which the above-described position 

matching is fully completed, the subtraction image Psu 
represents only the diseased portion K, which is not present 
in the first image Pi and is only actually present in the second 
image P2 and which is the actual difference between the images 

20 Pi and P2, as shown in the first image Pi (Fig. 6A) , the second 
image 2 (Fig. 6B) , and the subtraction image Psu (Fig. 6C) . 
However, it is difficult to obtain an actual perfect matching 
of the positions, so in actuality , not only the diseased portion, 
which is the actual difference between the images Pi and the 

25 image P2, appears in the subtraction image Psu, but also 
artifacts L caused by misalignment between the outlines of 
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structural elements such as bone tissue appearing in the two 
images Pi and P2 remain therein (refer to Fig. 7). 

At this point, the image processing means 12 of the 
anomalous shadow detecting apparatus 10 according to the current 
5 embodiment performs an image process that suppresses the 
artifacts 1 in relation to the diseased portion K on the 
subtraction image Psu. More specifically, first, the 
subtraction image Psu (Fig. 8A) is separated into a high-density 
image Psul formed of the portions of the image having a density 
O 10 equal to or higher than the base density (the entire image Psul 
j| is formed of said image portions having a density equal to or 

5 greater than the base density, and the low-density portions, 

St whi ch have been converted to the base densi ty) , and a 1 ow-densi ty 

J\ image Psu2 (Fig. 8B) formed of the portions of the subtraction 

jjj 15 image Psu having a density equal to or less than the base density 
« (the entire image PSU2 is formed of said image portions having 

H a density equal to or less than the base density, and the 

high-density portions, which have been converted to the based 
density). Continuing, the high-density image Psul is subjected 
20 to a morphology process (an opening process) employing a 
structuring element (round-shaped) of a predetermined size, 
(Fig. 8C). On the other hand, the low-density image Psu2 is 
subjected to a morphology process (a closing process) employing 
a structuring element (round-shaped) of a predetermined size. 
25 Here, as described above, the size of the structuring elements 
employed in the morphology process are of a size that is smaller 
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than the diseased portion K and larger than the width of the 
elongated artifacts L. The actual structuring elements to be 
employed are set based on the clinical data that has been gathered 
in advance for each case; in a case, for example, in which the 
5 size of a subtraction image Psu is 1780 X 1780 pixels, an element 
would be, for example, a round-shaped structuring element having 
a diameter of 4-5 pixels would be suitable. However, the size 
and shape of this structuring element can be changed according 
to the size and shape of the diseased portion that is to be 

O 10 extracted or the size and shape of the artifacts that are to 

g be eliminated. 

5§ Because the opening-process from among the aforementioned 

W several types of morphology processes is an operation consisting 

of eliminating the high-density image portions among the image 
15 portions smaller in size than the structuring element, aside 
2 from the diseased portion K, which is larger than the structuring 

^ element, the artifacts Ll (i.e., the artifacts having a 

high-density from among the artifacts L) that are smaller in 
size than the structuring element are substantially eliminated. 
20 in the same manner, the closing-process from among the 

aforementioned several types of morphology processes is an 
operation consisting of eliminating the low-density image 
portions among the image portions smaller in size than the 
structuring element, aside from the diseased portion K, which 
25 is larger than the structuring element, the artifacts L2 (i .e. , 
the artifacts having a low-density from among the artifacts L) 
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that are smaller in size than the structuring element are 
substantially eliminated (Fig. 8C) . 

The pixels of the subtraction image Psul and the 
subtraction image Psu2 obtained in this manner are again 
correlated, and both of said images are weighted at a 0.5:0.5 
ratio and subjected to an addition process to again form a 
processed combined subtraction image Psu' (Fig. 8D) . 

Because the processed subtraction image Psu' formed again 
in this manner is an image in which the artifacts L, which are 
caused by misalignments of the corresponding structural 
positions between the image Pi and the image P2 , have been 
suppressed (eliminated) relative to the diseased portion K, 
which is the actual difference between the image Pi and the image 
P2, the degree to which the artifacts appearing in a subtraction 
image interfere with the diagnostic reading thereof is reduced 
compared to currently available apparatuses. 

Because the artifacts 1 have been suppressed in this 
processed subtraction image Psu' , there is no hindrance to the 
diagnostic reading of the diseased portion K, and the diagnostic 
readability is improved in comparison to currently available 
methods and apparatuses. According to the current embodiment, 
the diseased portion K is again detected from this processed 
subtraction image psu' as an anomalous shadow by the detecting 
means 13. Hereinafter, the process performed by the detecting 
means 13 will be explained. 

First, the processed subtraction image Psu' is divided 
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into a plurality of small regions. At this time, the processed 
subtraction image Psu' is divided so that a portion of each small 
region AO overlaps, as shown in Fig. 9. The reason a portion 
of each small region AO has been caused to overlap is to prevent 
5 failure in the detecting of the anomalous shadow for cases in 
which the anomalous shadow, that is, the diseased portion is 
located at a border of any of said small regions AO. Note that 
in Fig. 9, in order to clearly show the state in which each small 
region AO overlaps, solid-line regions and broken-line regions 
O 10 are shown alternately. According to the current embodiment, the 
;U small regions AO have a size of 32 X 32 pixels. Here, in the 

2 early stages of lung cancer, the diseased portions thereof have 

[if a size in the 3-10 mm range. Because the size of the small regions 

!\ AO is 32 X 32 pixels according to the current embodiment, in 

2 15 order to detect the smallest diseased portion K, which has a 
Ji| size in the 3 mm range, the overlap of the small regions AO is 

|B * caused to be 3 mm (for a case in which the size of one pixel 

is 0.2 mm: 15 pixels) . 

Next, the number of pixels cn, which is the number of pixels 
20 having a pixel value above a threshold value Thl occurring within 
a small region AO, is counted, and the dispersion value a of 
the pixel value within a small region AO is added thereto. Then, 
for cases in which the ratio (%) of the number of pixels cn of 
the small region AO exceeds a predetermined threshold value Th2 
25 and the dispersion value a exceeds a threshold value Th3, that 
small region AO is selected as a small region Al in which there 
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is a high probability that an anomalous shadow is located 
therein. 

Here, if the pixel value is above the threshold value Thl, 
that pixel can be considered to represent an anomalous shadow; 
also, if said number of pixels cn is found to be above a certain 
percentage of the total number of pixels within the small region 
AO, there is a high probability that said pixels cn represents 
a suspected anomalous shadow. Further, for cases in which a small 
region AO is a so-called beta image, which has a uniform pixel 
value and appears not to contain an anomalous shadow, the 
dispersion value a of the pixel value becomes a small value 
because said small region AO contains only noise. Therefore, 
there is a high probability that a small region AO for which 
the dispersion value a does not reach the threshold value 3 is 
a beta image formed only of noise. Accordingly, a small region 
AO for which the number of pixels cn is above the threshold value 
Th2 and the dispersion value a is above the threshold value Th3 
has been selected as a small region Al of which there is a high 
probability that an anomalous shadow is contained therein. 

Here, because the processed subtraction image Psu' 
according to the current embodiment is an image in which the 
artifacts L have been suppressed, the accuracy with which a small 
region Al, which has a high probability of containing an 
anomalous shadow, can be detected is improved. 

in this way, when a small region Al has been selected, 
a determination is made as to whether or not a small round region 
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exists therein. This determination is made as follows. First, 
a center point c within the small region AO is obtained using 
only the pixels of which the pixel value is above the threshold 
value Thl. Note that Fig. 10 shows an image Gl formed of only 
5 the pixels within the small region AO whose pixel value exceeds 
the threshold value Thl. In this type of image, the center point 
c can be obtained as shown in Fig. 11. 

Next, taking the center point c as the center, the image 
Gl of the small region AO is converted to an extremes coordinate 
10 image. That is to say, the image Gl of the small region AO is 
converted to an extremes coordinate image P0 representing the 
distance r from the center point c and the angle 0 formed by 
the downward pointing vector passing through the center point 
c . 

15 Then, the edge portion extending in a horizontal direction, 

that is, in the direction of the angle 6 is detected from the 
extremes coordinate image P0, and an edge image E0 is formed. 
This edge image EO is formed by performing a convolution, by 
use of an edge mask M (Fig. 13) for detecting the horizontal 

20 edge portion, on the extremes coordinate image P0. An edge image 
EO is shown in Fig. 14. As shown in Fig. 14, the edge image EO 
has the pixel value of the pixels of edge portion El occurring 
within the extremes coordinate image P0, which is a pixel value 
of 0 (highest brightness) ; the pixel value of the other portions 

25 thereof have the largest pixel value (least brightness). 

Here, in order to improve the responsiveness of the edge 
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detecti ng mask m to the strai ght 1 i ne i n the hori zontal ch" recti on 
of the extremes coordinate image, the positive (+) portion of 
the mask is formed as an elongated elliptical shape so as to 
conform to the horizontal direction, and is selected so that 
5 the negative portion thereof is distributed on the upper side 
of the positive portion. In order to provide the mask with 
adequate orientational selectivity, a negative portion of this 
sort is indispensable. Then, the extremes coordinate image PO 
is subjected to a convolution using this mask, and an edge image 

10 EO is formed. 

when an edge image EO has been formed in this way, a 
histogram HO is formed in the vertical direction of the edge 
portion El occurring in the edge image EO, that is, in the 
direction r. As shown in Fig. 15, this histogram HO represents 

15 the distribution of the edge corresponding to the area from the 
center point c to the distance r. Therefore, as the edge is closer 
i n shape to a perfect ci rcl e , it becomes di stri buted equi di stant 
from the center point c, whereby the highest frequency of the 
histogram becomes large within said uniform distance. 

20 Accordingly, for cases in which the frequency of the histogram 
HO is higher than a threshold value Th4, it is determined that 
a substanti al 1 y round regi on i s present wi thi n the sel ected smal 1 
region Al. 

Next, the small substantially round region is detected. 
25 This detection is performed by matching subjecting the edge image 
EO to template matching by an elasticity template TO, which has 
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as its initial position the base template obtained by subjecting 
the round region formed of the concentration of points around 
a point 0, which is the highest frequency of the histogram HO 
at the radius r to extreme coordinate conversion. 

5 when this template matching is performed, all of the pixel s 

forming the elasticity template TO do not move as a unit; the 
template TO is a template that can receive the constraining of 
a virtual spring, wherein each pixel is constrained by the 
constraining power corresponding to the movement amount of each 

10 pixel in the interval between each said pixel and the pixels 
adjacent thereto, while being set so as to be capable of moving 
independently in the direction r, and the template is formed 
so that the entire template might be subjected to an elastic 
transform. Note that the base template is the straight line 

15 extending in the horizontal direction occurring at the highest 
frequency of the histogram HO. 

First, as shown in Figs. 16A and 16B, as to the initial 
position (the shape not having an elasticity transform) of the 
base template, the elasticity template TO is disposed above the 

20 edge image EO (Fig. 16a), and each of the pixels forming this 
elasticity template TO is made independent and moved in the 
direction r (the north-south direction occurring within the 
extremes coordinate image). At this time, the movement amount 
of each pixel of the elasticity template t' at this time is 

25 obtained as follows. First, for the peripheral range of each 
pixel (initial position ±r), the difference between each pixel 
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value p (nri+r, 0i) and the pixel value p at the initial position 
(nri, n6i). Here, the pixel value of the position at which the 
value of r is small is subtracted from the pixel value of the 
position at which the value of r is large. Then, the sum of this 
difference is obtained according to the formula (7) below: 
Sum difference of pixel values = 

{ g ( nri + r 3 n 6 i )- g ( nri ,n 6 i ) } ' 



N1 

20 

r 



N2 

+ 20 



{g(nn',n8i)-g(nri-r,n0i)} 
r 



(7) 



where Nl is a range having a larger r than the initial value 
among the peripheral ranges, and N2 is a range having a smaller 
r than the initial value among the peripheral ranges: 

x( x > 0) 



0 (x) = 



0 ( x < 0) 



This sum difference of the pixel values shows that a 
positive value is obtained if there are bright pixels (high pixel 
value, low density) in the direction having a value r larger 
than that of the initial position, and that a negative value 
is obtained if there if there are bright pixels (high pixel value, 
low density) in the direction having a value r smaller than that 
of the initial position. Further, by dividing the difference 
of the pixel values by r, the difference of the pixels in the 
vicinity of the initial position is weighted. That is to say, 
because the edge portion El is brighter than the periphery 
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thereof, if the edge portion El is located in a direction having 
an r larger than the initial position, a positive value is 
obtained for the sum difference of the pixel values, and if the 
edge portion El is located in a direction having an r smaller 
than the initial position, a negative value is obtained for the 
sum difference of the pixel values; this positive or negative 
symbol provi des an i ndi cati on of the ori entati on of the movement 
from the initial position, and an indication of the movement 
amount is provided by the absolute value thereof. 

At that point, the r of the movement amount (including 
the orientation) of each pixel of the elasticity template TO 
is defined by the formula (8) below, which utilizes a 
predetermined coefficient b: 

{g(nn' + r J n0i)-g(nri,n0i)} 



N1 
r 



N2 
r 



{g(nn\nfli)-g(nri-r,n0i)} ' 
r 



b 

•b (8) 



Although the movement amount rn obtained for each pixel 
20 in this manner, as shown in Fig. 16B, is the movement amount 
of each independently moved pixel , as described above, because 
this elasticity template TO is a template by which each pixel 
forming the image is constrained by the peripheral pixels of 
said elasticity template TO, the pixels are not moved the exact 
25 amount of the aforementioned movement amount rn, but are moved 
as determined according to each movement amount rk (k=n, n±l, 
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...) of the adjacent pixels (for example, both of the adjacent 
pixels: pixel n-1; pixel n+1), and also the pixels included up 
to said adjacent pixels (pixel n-2; pixel n-1; pixel n+1; pixel 
n+2; etc.), as shown in the Formula (9) below: 

5 

r n -2 ( a k" r k) (9) 

Here, it is preferable that the spring constant ak is set 
so as to be large with regard to the pixel of interest n and 
10 small with regard to the peripheral pixels n ± 1, .... That is 
to say, because the spring constant is defined as: 

a n >a n± i >a n± 2 > > a n±m (10) 

15 in formula (9), the movement amount rn of the pixel of interest 
n becomes the movement amount corresponding to the difference 
between the movement amount rn thereof and the movement amount 
rk of the pixels adjacent thereto, and is constrained by the 
virtual elasticity. 

20 By repeating the operation of moving each pixel of the 

elasticity template TO a little at a time as described above, 
the edge portion El can be detected with accuracy. Note that 
the determination as to when the repetition of the operation 
is to be terminated is based on the sum difference of the movement 

25 amount having fallen below a predetermined threshold value, or 
the number of repetitions having reached a predetermined number 
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of times . 

According to the process described above, a substantially 
round region is detected with a high degree of accuracy by the 
pixels forming the elasticity template TO, and finally, the 
5 pixels forming the elasticity template TO are returned to the 
actual image, the area between adjacent pixels is subjected to 
a correction process (a linear correction, a spline correction, 
etc.) , and by connecting the adjacent pixels by a closing curved 
line as shown in Fig. 17, the substantially round region can 
10 be extracted as a suspected anomalous shadow. 
Jit The image display means 20 displays the extracted 

X anomal ous shadow (the di seased porti on K) as an emphasi zed regi on 

in 

m enclosed in a closed region, as shown in Fig. 18A- Note that 

the anomalous shadow, as shown in Fig. 18B, can be not only 
Hi 15 enclosed within the closed region, but also can be specified 
S by an arrow indicator Y or the like. 

Because the artifacts have been suppressed and the 
diseased portion K enhanced in the image displayed on the image 
display means 20 in this way, the presence of the diseased portion 
20 can be more accurately recognized visually. 

According to the anomalous shadow detecting apparatus 
according to the current embodiment described in detail above, 
by subjecting an interimage image to an image process wherein 
the actual difference between the two images upon which said 
25 interimage is based is enhanced relative to the artifacts 

appearing therein, the artifacts, which have a negative effect 
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on the diagnostic reading of the interimage image, can be reduced 
relative to the actual difference between the two images upon 
which said interimage is based, which is the area concern therein 
Accordingly, the anomalous shadows can be detected with a high 
5 degree of accuracy with no negative effect due to the artifacts. 

Note that the according to the above-described embodiment, 
although an opening-process morphology process has been applied 
to the high-density subtraction image Psul, instead of this 
process, a highest-value filtering process can be applied, and 
l jJ 10 the same result obtained by application of said opening-process 
;J can be obtained thereby, in the same manner, although a 

Sj closing-process morphology process has been applied to the 

S low-density subtraction image Psu2, instead of this process, 

^ a lowest-value filtering process can be applied, and the same 

15 result obtained by application of said closing-process can be 
% obtained thereby. 

Further, because in the artifacts occurring in a 
concentration image the density changes radically in a narrow 
range, these sharply defined artifacts can be searched for by 
20 use of a Laplacian filtering, and by subjecting the artifacts 
obtained by said search to a smoothing process, said artifacts 
can be suppressed relative to the diseased portion, etc. actual 
difference between the two images on which said concentration 
image is based. In the same manner, because the artifacts include 
25 more high frequency components than the actual difference that 
represents a diseased portion or the like, the detected high 
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frequency components can be subjected to a blurring process. 
The blurring process to which the high-frequency components are 
subjected can consist of a Fourier transform followed by a 
deletion of the high-frequency components and a reverse Fourier 
5 transform, or a smoothing process employing a blurring mask. 
Note that for cases i n whi ch the above descri bed anomal ous shadow 
detecting apparatus according to the present invention is 
applied to a subject image such as a chest X-ray or the like, 
even if a blurring process is applied only in the vertical 

10 di rection, an adequate result according to the present invention 
can be obtained with respect to the artifacts caused by the bone 
structures extending mainly in the horizontal direction- 
Further, instead of performing a process to suppress the 
artifacts, the actual difference that represents a diseased 

15 portion or the like can be subjected to an enhancement process 
for enhancing said actual difference relative to the artifacts. 
For example, the low-frequency components forming the shadow 
of a tumor, etc. can be subjected to an enhancement process. 
Of course, it is permissible that both the suppressing of the 

20 artifacts and the enhancing of the actual difference be 
performed. 

Note that according to the current embodiment, upon being 
subjected to the non-linear transform (warping), as shown in 
Fig. 5, the first image Pi does not maintain the shape 
25 (rectangular) of the original image. Accordingly, the 

post-warping image Pi is not an image having a rectangular 
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outline such as that shown in Fig. 6A, but has an irregular 
outline formed of curved lines. On the other hand, because the 
second image P2 has not been subjected to any type of 
transformative process, said second image P2 has a rectangular 
5 outline, Accordingly, when correlating the pixels of both of 
said images and performing a subtraction process therebetween, 
there are cases for which there are pixels in the post-warping 
image Pi for which there are no corresponding pixels in the second 
image P2, and it is impossible to perform a subtraction process 

10 on such pixels. 

At that point, for cases in which there are pixels in the 
image that has been subjected to the transform process (in the 
current embodiment, the first image Pi) for which there are no 
corresponding pixels in the image that has not been subjected 

15 to the transform process (in the current embodiment, the second 
image P2) , said pixels for which corresponding pixels are not 
present can be reset and the value of said reset pixels can be 
caused to be the same as that of the pixels of the corresponding 
image (the second image P2) . 

20 Further, as a method instead of this pixel compensation 

method, the pixels of the image serving as the base image can 
be deleted at the same time as those of the image that has been 
subjected to the transform process. This is due the fact that 
there are extremely few cases in which the image portions for 

25 which there are no corresponding pixels represent an area of 
interest with respect to the diagnostic reading, and said pixels 
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can therefore be considered as unnecessary. Further developing 
this approach, only the area of interest in each of respective 
post-warping image (the first image Pi) and bas image (the second 
image P2) can be extracted, so as to limit in advance the area 
5 of the images to be subjected to the subtraction process, in 
this manner, the computation amount is advantageously reduced 
by limiting in advance the area of the images to be subjected 
to the subtraction process. According to the embodiment 
described above, because chest X-ray images have been applied, 

10 it is preferable that a thorax recognition process, whereby all 
of the pixels representing portions outside of the chest cavity 
are assigned a uniform value, is performed on the post-warping 
image (the first image Pi) and the base image (the second image 
P2) . Note that, as to the thorax recognition process, the 

15 technology disclosed in Japanese Unexamined Patent Publication 
No. 8(1996)-335271, etc. can be applied. 

Still further, according to the above-described 
embodiment, although the detecting means 13 divides the 
processed subtraction image Psu 1 into small regions having a 

20 size of 32 X 32 pixels, anomalous shadows occur in a variety 
of sizes; therefore, the processed subtraction image Psu' can 
be divided into series of small regions, each of a different 
size, for each iteration of which the above described process 
for detecting the small round regions occurring therein is 

25 performed, whereby various sizes of anomalous shadows can be 
detected . 
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in addition, according to the above-described embodiment, 
although substantially round shaped diseased portion K has been 
detected as an anomalous shadow, by use of a substantially 
rectangular shaped structuring element instead of a 
substantially round structural element by the image processing 
means 12 when obtaining the processed subtraction image Psu', 
the diseased portions having a substantially rectangular shape 
can be compared to the artifacts and enhanced relatively thereto, 
in this case, by the subjecting of the processed subtraction 
image Psu' to a template matching process employing a 
substantially rectangular template and performed by the 
detecting means 13, the substantially rectangular anomalous 
shadows can be detected. 
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